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AFIT/GE/ENG/93D-31 


Abstract 

The resolution of optical imaging systems is severely degraded from the diffraction limit by the 
random effects of the atmosphere. Techniques exist to compensate for the atmospheric turbulence, 
one of which is adaptive optics. A critical component in the adaptive optics system is the wavefront 
sensor. Presently, two types of sensors are being used-the Hartmann-Shack Wavefront Sensor and 
the Shearing Interferometer. Previous studies have compared these two sensors and found them to 
perform identically for a point source. However, to date, no comparison has been performed for an 
extended source and subaperture spacing larger than the correlation length of the atmosphere, r a . 
This thesis has examined this problem and compared these two sensors for the above conditions. 
Results indicate that the sensors have comparable performance when an infinite number of photons 
are available. However, the photon limited cases indicate superior performance by the Hartmann 
Sensor. 
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PERFORMANCE COMPARISON OF SHEARING INTERFEROMETER 


AND 

HARTMANN WAVE FRONT SENSORS 

I. Introduction 

The upper bound on the resolution achievable by a ground based imaging system is set by the 
diffraction limit. This diffraction limit is inversely proportional to the 9ize of the optical system. 
However, Sir Isaac Newton recognized in 1730 that another, more restrictive limit existed when he 
wrote: 


If the Theory of making Telescopes could at length be fully brought into Practice, 
yet there would be certain Bounds beyond which Telescopes could not perform. For the 
Air through which we look upon the Stars, is in perpetual Tremor; as may be seen by 
the tremulous Motion of Shadows cast from high Towers, and by the twinkling of the 
fix’d Stars [16]. 

This “Tremor” of the air is today known as atmospheric distortion and has been the object of 
considerable research for the past several decades [2], [9], [10], [14], [19]. 

1.1 Atmospheric Turbulence and Its Effect 

Atmospheric turbulence arises from the sun’s uneven heating of the Earth. This differential 
heating generates pockets of air, called eddies, each a different temperature. Each eddy deflects, 
or refracts, the light by an amount related to its temperature. Hence, as light passes from one 
air pocket to the next, it is randomly bent as the temperature changes. As the entire wavefront 
propagates through the atmosphere, it passes through different “eddies” and hence is bent by 
different amounts, (see Figure 1.1), creating uneven lines of constant phase. The resulting image 
formed by a telescope looking through the atmosphere is distorted and blurred. 
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Figure 1.1. Light Experiencing Atmospheric Turbulence. 


From a communication theory perspective, the atmospheric turbulence can be thought of 
as a low pass filter with a spatial frequency cut-off far below the diffraction limit. And just as 
sharp transitions, such as the corners of a square wave, require high frequencies to pass through 
a communication channel, so do the sharp transitions that might be present on a satellite in 
Earth orbit. The atmospheric low pass filter effectively removes these high frequency components, 
resulting in a blurred image. This cut-off frequency is related to Fried’s Coherence Diameter, r 0 [2], 
which can be thought of as the maximum separation of two points for which the field is correlated. 
This parameter is, of course, related to the strength of the atmospheric turbulence. 

1.2 Adaptive Optic Compensation 

In 1953, Horace W. Babcock [1] proposed using a mirror coated with an oil film to correct 
for the phase distortions imposed by atmospheric turbulence. The thickness of the oil film would 
be altered by electric charge, changing the optical path length such that the effects of atmospheric 
turbulence would be canceled out. This idea is the basis of current adaptive optic systems. Present 
day systems are composed of three major components [15]: a phase-shifting optical element, a 
wavefront sensor, and a servo-control system, (see Figure 1.2). As the light passes through the 
telescope, the atmospheric phase distortions are compensated for by using the phase-shifting optical 






Figure 1.2. Block Diagram of a General Adaptive Optic System. 

element, which is usually a deformable mirror. The reflected light is then analyzed with the 
wavefront sensor which detects the phase along the wavefront. This information is used by the 
servo-control system, usually a dedicated computer or specialized hardware, to compute the residual 
phase error and control the deformable mirror, thus closing the loop. Systems using this general 
approach have been under extensive development for the past twenty to thirty years. Recent 
declassification of this technology has revealed impressive and tantalizing results [3], [4], [17]. 

1.3 Measurement of Wavefront Phase 

The phase of the incident optical wave cannot be measured directly. One method of indirectly 
measuring the wavefront phase is to measure the phase slope of small zones of the wavefront. These 
zones, within the overall aperture of the telescope, will be referred to as “subapertures.” These 
measurements can be used to create a “map” of the slope and hence the phase along the entire 
wavefront. An illustration of how such a phase map might be constructed is given in Figure 1.3 [19]. 
Here, s x or s v is the phase slope in the x or y direction, respectively, obtained from the wavefront 
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Figure 1.3. Phase Map Reconstruction. 

sensor. 4>ij is the wavefront phase at location i,j. Therefore, the phase slope could be written as: 

x _ ^0,1 - <*-1,1 

51 - - d 

y _ — 0 - 1,0 

Sl -- 

( 1 - 1 ) 

Examination of Figure 1.3 shows that there are more slope measurements than nodes where we 
wish to determine the phase. Thus, a least squares algorithm can be used to reconstruct the phase 
wavefront map [10]. 

Two instruments are commonly used to measure the phase slopes: the Hartmann-Shack 
Wavefront sensor and the Shearing Interferometer Wavefront sensor. The Hartmann-Shack Wave- 
front sensor, (see Figure 1.4), measures the slope by focusing the wavefront incident onto each 
subaperture to form a spot. The slope of the wavefront segment is proportional to the offset of the 
spot from center. The Shearing Interferometer, (see Figure 1.5), determines the wavefront slope 
by duplicating the wavefront and shifting, or shearing, it with respect to the original. The two 
wavefronts are interfered to produce an interferogram. The slope information along the wavefront 
can be determined from the information in the interferogram. The Hartmann Wavefront sensor 
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X shear 


Figure 1.5. Diagram of the Shearing Interferometer. 
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and Shearing Interferometer, and how each one measures the phase slopes, are discussed in greater 
detail in Chapter II. 

As with all practical systems, the Hartman Wavefront sensor and Shearing Interferometer do 
not measure the wavefront perfectly. Two significant, and unavoidable, errors associated with wave- 
front sensing are photon noise effects and sampling errors. As the brightness of the source decreases, 
photon noise can become a considerable problem. Photon noise results from the uncertainty in the 
detection of the light due to fluctuations in the photo count statistic, often modeled as a Poisson 
process [8:pg. 482]. Photon noise can be decreased by either making the light source brighter, 
increasing the integration time for each sample or making the subaperture larger. Sampling errors 
occur due to the finite size of the subaperture. Instead of measuring the exact wavefront slope 
within each subaperture, the average wavefront slope, over the subaperture, is measured. Sampling 
errors can be reduced by either making the subapertures smaller or decreasing the integration time 
for each sample. A tradeoff must therefore be made between photon noise effects and sampling 
errors, the optimum condition being when these errors are equal [10:pg. 687]. 

Previous studies have shown that both of these sensors perform equally well under quasi-ideal 
conditions, i.e. when imaging a point source and there is little or no atmospheric turbulence [10]. 
However, no studies have been performed comparing the sensors when using an extended source, 
and the atmospheric turbulence is such that the correlation length, r„, is less than the subaperture 
diameter. As the source becomes larger and r 0 smaller, the randomness of the wavefront phase 
increases-to a point where the phase at one end of the subaperture is not reiated to the phase 
at the other end of the subaperture. As r„ decreases, the phase will vary within the subaperture 
diameter. Thus, the sensor does not “see” the true tilt, (called z-tilt), but instead “sees” the 
average of the true tilt over the subaperture diameter, (called g or s tilt). Phase measurements 
under these conditions will not reflect the true wavefront tilt. 
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This thesis will examine the sensors under these conditions and develop mathematical models 
to predict their performance. 

1-4 Scope 

The scope of this thesis will be to develop an expression to model the performance of each 
of these sensors given the subaperture size, light level, source size and correlation length, r 0 . The 
performance of the two sensor will be compared based upon these models. 

1.5 Summary of Key Results 

This analysis determines that, as expected, the Hartmann Wavefront Sensor and Shearing 
Interferometer perform comparably for weak atmospheric turbulence, (large r n ), and very small 
sources, or very bright sources. However, for strong atmospheric turbulence, extended sources and 
photon limited conditions, the Hartmann Wavefront Sensor displays better performance. 

1.6 Overview 

Chapter II briefly summarizes the theory involved with wavefront sensing. Chapter III ex¬ 
amines the Hartmann-Shack Wavefront sensor. Performance models for the Hartmann WFS are 
developed. Chapter IV discusses the Shearing Interferometer Wavefront sensor and performance 
models are developed. Chapter V compares the sensors based upon the models developed in Chap¬ 
ters III and IV. Results are presented that indicate which sensor performs the best under various 
situations. Chapter VI gives some conclusions about the results and recommendations for further 
study. 
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II. Preliminaries 


Before proceeding with an analysis of the errors resulting from wavefront measurements with 
the Hartmann Wavefront sensor and the Shearing Interferometer, it will be useful to establish 
how each sensor measures the wavefront. The Hartmann Wavefront sensor will be discussed in 
Section 2.1 and the Shearing Interferometer in Section 2.2. Finally, Section 2.3 establishes some 
definitions that will be used throughout this thesis. 


2.1 Hartmann Wavefront Sensor 

A Hartmann Wavefront sensor is composed of an array of subapertures which segment the 
wavefront. In each subaperture is a lens which focuses the segmented wavefront onto an array of 
photo detectors that determine the resulting spot positions. A single element of the lens array, 
illustrated in Figure 1.4, is depicted in Figure 2.1. By use of simple geometry, we can see that the 
focused spot centroid, z e , is easily related to the tilt, a, of the incoming optical wave: 


tan(a) = —■ 
ft 



( 2 . 1 ) 


However, the question as to how the centroid of the intensity pattern is determined still 
remains. Using a detector array like the one in Figure 2.2, Kane [12] determined that the centroid 


Plane wave 
with no bit 




Figure 2.1. One Subaperture of the Hartmann Sensor with (a) a Normally Incident Plane Wave 
and (b) a Tilted Plane Wave. 
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Figure 2.2. Model of Photo Detectors in Focal Plane of the Hartmann Sensor. 




(a) (b) 

Figure 2.3. (a) Array of Spots from a Plane Wave, (b) Array of Spots from a Distorted Wave. 

of the spot could be estimated using: 


U{xt,y e ) = 


XjPjj 

StE;Pij 


( 2 . 2 ) 


where U(x e , y e ) is the estimate of the centroid in the x direction, x,- is the center of the x* h detector 
pixel and j is the number of photon events in the (*, j) th pixel. Note that Equation 2.2 looks very 
much like a one dimensional centroid definition. Figures 2.1 and 2.2 represented one subaperture. 
If an array of subapertures is considered, a pattern much like Figure 2.3 might be seen where each 
spot results from one subaperture. Figure 2.3(a) represents the pattern that might be seen in the 
detector plane of a Hartmann Sensor illuminated by a plane wave with no atmospheric distortion. 
Figure 2.3(b) represents the same pattern resulting from a distorted wavefront. 
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2.2 Shearing Interferometer 


As discussed in Chapter I, the Shearing Interferometer takes a different approach to wavefront 
sensing. Instead of spatially segmenting the wavefront, the Shearing Interferometer duplicates 
the incident wavefront; shifts (shears) it; and interferes it with itself. Several types of Shearing 
Interferometers exist [11]. For this thesis, a version of the AC Shearing Interferometer is considered- 
which induces a temporal modulation on one of the sheared beams by use of a rotating diffraction 
grating [ll:pg- 183]. The temporal modulation causes the interference fringes to evolve with time in 
the subaperture. The phase of the fringes can then be determined from the temporally modulated 
wavefront. Exactly how the interference fringes are related to the wavefront distortion is discussed 
in greater detail in Section 2.2.1. Section 2.2.2 describes how the phase can be derived from this 
intensity pattern using Goodman’s Discrete Fourier Transform (DFT) Fringe Estimator [8:pp. 494- 
501]. 


2.2.1 Relationship Between Detected Intensity and Wavefront Tilt for the Shearing Inter¬ 
ferometer Imperative to understanding how the Shearing Interferometer works is determining how 
the interference pattern is related to the wavefront tilt. The incident wave, before being sheared, 
can be expressed as: 

U,(x,y) = E(x,y)e>«*-*> (2.3) 

where E(x, y) is the field from the source and <f{x, y) is the phase distortion due to atmospheric 
turbulence. The field from the source, E, is assumed to be complex, so that <f>(x, y) is only the 
phase due to atmospheric distortion. Recall that the resulting intensity pattern will be sampled 
temporally. Thus, it is reasonable to assume that the integration time for each sample will be short 
enough such that Taylor’s frozen turbulence hypothesis can be assumed [8:pg. 433]. Hence, <j>(x, y) 
will be assumed to be deterministic. The two sheared wavefronts are written as: 

Ui(z,y) = Uj ^x + ^r,y) and U 2 (x,y,<) = U, (x - ^-.v) * >wt (2-4) 
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where the exp(juf) designates that a temporal modulation has been applied to U 2 using a rotating 
grating. The total field in the detector plane is: 


U, = Ui+U 2 (2.5) 

The intensity pattern in the detector plane of the Shearing Interferometer can be described as: 

/«(*,{/,<) = <U«U;) (2.6) 

where the a* indicates the complex conjugate of a and the t dependence results from the ac mod¬ 
ulation applied to U 2 . Substituting Equation 2.5 into Equation 2.6 and expanding gives: 

/.(*, y, t) = (|Ui| 2 ) + (|u 2 | 2 ) + (UiU; + u;u 2 ) (2.7) 

The first two terms are simply: 

(|U X | 2 ) = (|U 2 | 2 ) = (|£(*,y)| 2 ) = 7(*,y) (2.8) 

Therefore, substituting Equation 2.8 into Equation 2.7 we obtain: 

/,(*, y, t) = 2I(x, y) + (UiU; + UJU 2 ) (2.9) 

Now looking at the cross terms, 

(U X U; + U X U 2 ) = ^E (x + y) (x - y-, y) *■>)->»* 

+E* (*+ ^.y) e-^( ,+ ^-*)E (x- ^.y) (2.10) 
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Which simplifies to: 


(UiUS + UtU,) = ^E (*+ E* (x - ^,y) 

+£* (x+ E ^x - ^,y) e*W—'¥.»)-♦(*+¥•»)+««]^ ( 2 .ll) 

Note that since <t>(x, y) is deterministic, the expected value operators can be pulled in to include 
only the source field terms. 

(UiU? + UtU 2 ) = ^E (x + ^,y) E* (x - ^,y^ ■*)+««] 

+ ^E* (x + Y". v) E (x - y) ^ (2.12) 

To simplify the notation, define: 

A<t>(x, y, Ax) = 4> (z - y^ - 4> (z + y^ (2.13) 

Further, referencing Goodman [8:pg. 174], the mutual coherence function can be defined for this 
case as: 

T(Ax,0) = (e- (x + ^,y) E (x - y)^ (2.14) 

Thus Equation 2.12 becomes: 

(Uxu; + u;u 2 ) = r(Ax, o) e & A «*'»' A *> +< *'*i + r*(Ax, o)«- iIA ** > * iA *> + “ , i (215) 

Again referencing Goodman [8:pg. 180], for quasimonochromatic conditions, the mutual coherence 
function can be written as: 

T(Ax, Ay) = J(Ax, Ay) (2.16) 
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where J is the mutual intensity and can be found from the Van Cittert-Zernike Theorem. From 
Van Cittert Zernike, we can relate J to the intensity distribution of the source, /,, by [8:pg. 209]: 

J(Ax, Ay) = K JJ exp 

where K is some, possibly complex, constant. Thus, J can be found by taking the Fourier Transform 
of the intensity distribution of the source. It should be noted that a Shearing Interferometer will 
consist of two sections. One will shear in the Az direction and the other in the Ay direction. As 
both legs are otherwise identical, only the Az shear will be considered and Ay can be assumed to 
be zero. Further note that for the situations of interest, the source distributions will be assumed 
to be symmetric about the center of the source, (e.g. circular, rectangular or Gaussian). Thus by 
Fourier Transform properties [6:Table 7-1], if /(z) and F(C) are Fourier Transform pairs, and if 
f(x) is real and even, then F(Q will be real and even. Hence, for the situations of interest, J is 
assumed to be real and 

J* = J = / (2.18) 

Thus, Equation 2.15 can be written as: 


2x 

J=~ (AxC + Ayr;) 
A z 


dCd»? 


(2.17) 


(UlU; + u;u 2 ) = J(Ax,0) [«bA**.».A*)+wtl + e -;[A*«,v,A*)+<w]j 

= 7(Az,0)[2cos(A^(z,y, Az) + urt)] (2.19) 

Here, I will continue to follow Goodman and define the normalized mutual coherence function, or 
degree of coherence function, as: 


/A A x r(Az, Ay) 


( 2 . 20 ) 
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Recall that: 


7(0,0) = ^|£(z,y)| 2 ) = /(r,v) 


( 2 . 21 ) 


and Equation 2.20 can be written as: 

J(Ax, 0) = /(*, y)fi(Ax, 0) (2.22) 

Substituting Equation 2.22 into Equation 2.19 yields: 

(UiUJ + UJU,) = 2/(*,y)/»(A*,0)cos(A^(*,y,A*) + «0 (2.23) 

If the fields Ui and U 2 are assumed to have equal amplitude, then the complex coherence factor 
is equal to the visibility, V, of the fringes formed by the interference [8:pg. 163]. Equation 2.23 can 
be substituted back into Equation 2.9 to yield an intensity pattern of: 

I*(x, y, t ) = 2 1(x, y) [1 + V cos ( A+{x, y, Ax) + wt)] (2.24) 

Equation 2.24 shows that the wavefront difference, A<p, is encoded in the phase of the ac modulated 
fringe pattern. Therefore, to determine the wavefront tilt, the phase of the fringe pattern needs to 
be determined. That process is discussed in Section 2.2.2. 

2.2.2 Shearing Interferometer Detector Model The problem of retrieving the phase from an 
ac signal, such as Equation 2.24, is a common one in communication systems [18]. One method is 
the Discrete Fourier Transform (DFT) Estimator [20], [8:pp. 494-501]. Before describing the DFT 
Estimator in detail, a discrete representation of the detector plane intensity pattern, Equation 2.24, 
is needed. Equation 2.24 can be written as: 

I t (x,y,t) = 2I 1 + Vcos + A^(x,y,Air)^| (2.25) 
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where T is the temporal period of the fringe pattern and I is the intensity of each sheared beam in 
the detector plane. The detector plane is segmented, (as for the Hartmann sensor), into discrete 
detectors or subapertures. Here, it is assumed that each subaperture is small enough such that 
/(x, y) is essentially constant across the subaperture. The intensity in each of these subapertures, 
defined by W,(x,y), will be integrated for a time, r, to form a vector of discrete photon count 
values, see Figure 2.4. Through this count vector, and Goodman’s DFT Fringe Estimator, the 
phase of the fringes can be determined. 

Using Equation 2.25, the count vector, K(n), could be described as: 


K (n) = a JjJ I(x, y, Ax, t)W,(x, y)W t (f)dfdxdy 


(2.26) 


where a is a proportionality constant that converts energy to photo events and is given by [8:Equation 
9.1-9], W, is the aperture function defining the detector pixel and W t is the temporal windowing 
function of each measurement interval. W,(x,y) and W t can be described by: 


W,(x,y) = rect rect 0) = 


W t (t) = rect 


(-)-{: 


- § < X < I and - f < y < f 
otherwise 

rn~^<x<rn+^ 

otherwise 


(2.27) 


(2.28) 


where d is the width of the detector and r is the integration time for each intensity measurement, 
n. Substituting Equation 2.25 into Equation 2.26 yields: 

K(n) = 2 la + V cos + A0(x, y, Ax)j J W,(x, y)lV)(t)dfdxdy (2.29) 
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I<x) 



Figure 2.4. Fringe Pattern being Sampled Temporally. 


Equation 2.29 is evaluated in Appendix B and shown to be: 

K(n) = 2 IarA + V cos + 2A*(A*) j j (2.30) 

where A<£ was defined as: 

AA^(Az) = Jl A^(jf,y,Ax)W,(r,y)dxdy (2-31) 

and the relationship: 

r-N = p 0 T (2.32) 

was used where r is the integration time, N is the number of samples in the count vector, K(n), i.e. 
0 < n < N — 1, p« is the number of fringes embraced by the count vector and T is the period of the 
fringe. Equation 2.30 describes the photon count vector. This is, of course, a discrete representation 
of the fringe pattern. However, we still need to determine A<f> from this count vector. 

The Shearing Interferometer can be set up such that p 0 is an integer number of fringes. The 
phase difference, A<£, can then be determined from K(n) using Goodman’s DFT fringe estimator. 
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The DFT of Equation 2.30 is: 


1 N ~ l 

H(p) = ^ Y1 ^(n)e ,(2 '" p/A,) (2.33) 

ASO 

Define the real and imaginary components of the DFT evaluated at p, as: 

Hr = fte{W(p.)} = ^ E *(») coe (2-34) 

nsO 

Hi = Im{7f(p 0 )} = i £tf(n)sin^ (2.35) 

ASO 

Using Equation 2.30, the mean values of Equations 2.34 and 2.35 can be evaluated to be: 


Hr = 

2IaAr - x 

2 Vcos(A^) 

(2.36) 

Hj = 

2IaAT ,— 


—2—Vsin(A^) 

(2.37) 

From which we can easily see that: 



= tan(A^) 

Hr 


(2.38) 

It should be noted that my model is slightly different from Goodman’s. 

I have assumed the 


fringes are modulated and hence pass by a single detector. Goodman assumed that the fringes are 
stationary and an array of detectors spatially sample the fringe pattern. However, Equation 2.30 
matches Goodman’s Equation (9.4-2). Therefore, the statistics Goodman develops for his DFT 
Estimator are valid for my model as well. These statistics will be used in Chapter IV. 
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2.3 Definitions 


2.3.1 z-tilt z-tilt, & t , is the wavefront slope, or tilt, over a subaperture, z-tilt is defined 
as [24]: 

S iK-g • * + y ■ y) y)W{x, y)dxdy 

^kJ[(x 2 + y 3 )lV(x,y)dxdy (239) 

where W(x,y) defines the subaperture over which we are considering the wavefront tilt; x and y 
are unit vectors in orthogonal directions; k = 2r/A, where A is the optical wavelength; and <f> is the 
turbulence induced phase distortion. 


2.3.2 g-tilt 3 -tilt, 9 e , is the wavefront slope, or tilt, as measured by the Hartmann-Shack 
Wavefront sensor. It is directly related to the instantaneous angular position of the centroid of the 
intensity distribution in the focal plane of the Hartmann Sensor. 3 -tilt is defined as [24]: 


0 y ) \ w (*> y)l 2 dxd y 

ijfllWXz, y )| 2 dxdy 


(2.40) 


where V is the gradient operator. 

2.3.3 s-tilt s-tilt, 9,, is the wavefront slope, or tilt, as measured by the Shearing Interfer¬ 
ometer. 9 S is defined by: 


9 , = 


y)dxdy + yff *«‘^ W(x, y)dxdy 


kjfW{x,y)dxdy 


(2.41) 


where Ax and Ay is the amount of shear in the x and y direction respectively. Note, in Equa- 
tion 2.41, that as Ax, Ay —► 0 then 9, —► 9 e defined in Equation 2.40. 
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2.3.4 Subaperture The Hartmann Wavefront sensor and the Shearing Interferometer mea¬ 


sures the wavefront tilt over some finite area. That area is the subaperture which is defined as: 





W(x, y) = rect 

;) (; 

i)- 


— | < x < | and 
otherwise 


!<»<! 


(2-42) 


With this subaperture definition, several items about Equations 2.39, 2.40 and 2.41 can be noted. 
First, in Equation 2.40, since fV(x, y) is real and only takes values of one or zero: 


m*,y)F = W{x,y) 


(2.43) 


Further the denominators of Equations 2.39, 2.40 and 2.41 can be evaluated as: 

\ k jj (* 2 + y 2 ) W{x,y)dxdy = 
k jj W(x, y)dxdy = kd 2 

and 

k jj W(x,y)dxdy = kd 2 , 

respectively. 


(2.44) 


(2.45) 


(2.46) 


2.3.5 Normalized Shear In the Shearing Interferometer, the wavefront is sheared by a dis¬ 
tance Ax before being interfered with itself. Let the normalized shear be defined as: 


s = 


Ax 

r 0 


(2.47) 


where r 0 is the atmospheric correlation length. 
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III. Hartmann-Shack Wavefront Sensor 


Two major sources of error have been assumed for the Hartmann Wavefront Sensor. Those 
two sources can be described by the error between and z-tilt in the wavefront measurement 
and the photon or shot noise effects. These two effects are independent and hence can simply be 
added to give the total mean square (MS) error. For notational simplicity, the total MS error for 
the Hartmann sensor will be defined as: 

Eh = Etfp + Eh } , (3.1) 

where Eh p is the MS error resulting from the detectors inability to accurately determine the true 
centroid due to photon noise effects and Eh 3 z is the MS error due to the difference between the 
true wavefront tilt, z-tilt, and the tilt seen by the Hartmann Sensor 7 -tilt. Eh will be manipulated 
to have units of (radians/r D cell) 2 . 

3.1 Shot Noise, Eh p 

An optical wave passing through a positive lens will be Fourier transformed at the lens’ back 
focal plane. Thus, one would assume that an incident plane wave will form an impulse function at 
the back focal plane. However, this would require the lens to have an infinite spatial bandwidth, 
which is not the case. Hence we must take diffraction effects into account. Diffraction will result 
in a plane wave being transformed to a finite width spot. Kane [12] has shown that the width of 
this spot is related to the MS error in determining the spot centroid by: 



where x r . is the true centroid location, x e is the centroid estimate given by Equation 2.2, rj is related 
to the detector efficiency, N is the total photon count and <r,- is the RMS width of the intensity 
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pattern in the detector plane. <r,- will be determined for two cases. Section 3.1.1 will determine the 
RMS spot width for a point source and Section 3.1.2 will determine for am extended source. 

3.1.1 RMS Spot Width for a Potnt Source The RMS spot width for a point source is needed 
for two reasons. First, as indicated by Equation 3.2, the MS error in determining the spot centroid 
is proportional to the RMS spot width. By finding the RMS width for a point source, the MS shot 
noise error, due solely to diffraction, can be determined. Second, the intensity pattern resulting 
from a point source also gives the diffraction limited Point Spread Function, (PSF), for a Hartmann 
subaperture. This relationship will be used later in Section 3.1.2. 

The diffraction pattern for a square lens of width d is given by [7: pg. 63]: 

«*■*»- <33) 

where x a and y„ are points in the focal plane and sinc(x) is defined as sinc(x) = sin(jrx)/xx [6:pg. 
45]. By defining: 

dx 0 , dy 0 

“ITT “ d (34) 

and normalizing I to one at x, y = 0, Equation 3.3 can be written as: 

I'(x,y) = sinc 2 (x)sinc 2 (y) (3.5) 

which represents the normalized intensity pattern for a point source with no atmospheric distor¬ 
tion. Analytical and numerical computations with a sine 2 function can be difficult. To ease the 
computations, the approach taken by Kane [12] is followed, and this distribution is approximated 
by a normalized Gaussian function having the same width at the e -1 points. 

Examination of Figure 3.1(a) shows that I'(x) has a value of e _1 at x = 0.5234. A normalized 
Gaussian curve can be matched to the sine 2 pattern by finding a variance, <r 2 , that forces the 
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normalized Gaussian function to have a value of e _1 at x = 0.5234. The Gaussian curve is given 
by: 

(3 ' 6) 

and the normalized Gaussian function by: 


G'(x,ir) 


G(x, tr) 
G(0,<r) 


= exp 



(3.7) 


It can easily be shown that a variance of a 2 = (0.37) 2 will match the two fimctions at the e~ l points. 
A normalized Gaussian function with such a variance is plotted in Figure 3.1(b) for comparison. 
Accounting for the variable change in Equation 3.4, the approximated intensity pattern can be 



which also represents the diffraction limited PSF for the Hartmann subaperture. 


To determine the shot noise error, En p , the RMS width of the approximated intensity dis¬ 
tribution, described by Equation 3.8, needs to be determined. However, since Equation 3.8 is a 
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Gaussian function, the RMS width is simply its standard deviation; hence: 


<Ti = 0.37^- 
a 


This result is similar to one obtained by Kane [13] for a circular aperture. 


3.1.2 RMS spot width due to Extended Source The analysis in Section 3.1.1 gave the RMS 
spot width for a point source. However, we are also interested in the effects due to an extended 
source. The source intensity, I ,, is modeled as an elliptical Gaussian distribution. This distribution 
is expected to be formed from a laser guide star [5:pg 1722]. 


T*(C,»?) = hex p 


L<i±M' 

L 2 «r? 


(3.10) 


where the constant b is added in to make the distribution “elliptical.” The intensity in the de¬ 
tector plane is given by the convolution of the source intensity with the diffraction limited PSF, 
(Equation 3.8): 


r / ' j \ x ' 2 + t / 2 x ' 3 + ( V ) 2 

Id(x , V) = exp- 2 * exP-—— r 


2(0.37)2 (*£)* 


M 


(3.11) 


where the constants have been ignored. This is essentially the convolution of two 2-D Gaussian 
functions, which can easily be evaluated to be: 


h{x\ ]/) = exp - 


& 

2 (4 £ +( z ‘f t ) 2 ) 2 ^+(if *.yy 


(3.12) 


where: 


= (0-37) 2 


(3.13) 
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is the diffraction limited spot width for a square aperture, determined in Section 3.1.1, and the 
factor of fijz converts the source RMS width, a,, from object space to image space. Thus we see 
that the intensity distribution in the detector plane of the Hartmann sensor is Gaussian with MS 
width equal to the sum of the diffraction limited MS width and the MS width of the source in 
image space. Hence the RMS spot width is given by: 

<Ti, = \J*dl + = (3-14) 


and in general: 

= (0 37) 2 (~p) 2 + (jfj (315) 

3.1.3 Mean Square Error Due to Finite Spot Size Recalling Equation 3.2, the RMS spot 
width, given by Equation 3.14, is related to a centroid location MS error. Using simple geome¬ 
try, (see Figure 3.2), this location MS error can be related to a MS error in the wavefront slope 
measurement by: 



(3.16) 


where is given by Equation 3.2, A <t> is assumed to represent the change in phase of the wavefront 
over the subaperture and the factor of two accounts for an identical error in the y direction. Note 
that <n has units of meters, as does fi and d (or meters per diameter). Thus Equation 3.16 has units 
of (meter/diameter) 2 . To transform these units to (radians/r 0 cell) 2 , a multiplication by ( kr 0 /d ) 2 , 
which has uu'ts 



( radians^ 2 ( meter \ 2 ( diameter^ 2 

( radians • diameter^ 2 

V meter ) \r 0 cell,/ \ meter J 

\ meter • r 0 cell ) 


(3.17) 
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Figure 3.2. Diagram of a Hartmann Sensor Subaperture 


is required. Thus 



and the shot noise MS error for the Hartmann Wavefront Sensor can be defined as: 


(3.18) 



(3.19) 


Substituting in Equation 3.15 for of and simplifying, we obtain: 



(3.20) 


For analysis of arbitrary subaperture sizes and seeing conditions, it is important to note that the 
second term in Equation 3.20 is only valid for d/r 0 < 1, where r 0 is Fried’s coherence diameter [2]. 
The term 0.37/jA/d represents the diffraction limited spot size in the detector plane. However, 
as the diameter, d , gets larger, the spot size does not shrink in a proportional way—instead, it is 
limited by the atmospheric correlation length, r„ [8:pp. 431-432]. Thus, for d/r 0 > 1, the d in the 
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second term needs to be replaced by a r a . Equation 3.20 can more accurately be written as: 


j *(*)' 

(¥)’<’?+(^)’ 

w 

(to)’^ + (0.37/,A) J 


for < 1 

for -£ > 1 
“• 


Performing algebraic simplification, recalling that k = 2t/A, Equation 3.21 becomes: 


(3.21) 


(3.22) 

\ 2 V[(V ^) J + (0.37) 2 ] for^>l 

Note that <r,/z is equivalent to the angular size of the source, while A/r, is the angular resolution 
limit induced by the atmosphere. The ratio of the angular size of the source to the angular resolution 
of the unaided optical system is defined by: 


£l 


t- z A 


(3.23) 


Using Equation 3.23, Equation 3.22 can be written as: 


Eap = { 


2 !»£[^ + (0.37)’(!*)’) 
( 2^C[/3> + (0.37)’] 


for f- < 1 

r a 


for A > 1 
~0 


(3.24) 


3.2 Mean Square Error between z-titt and g-tiH, Eg t , 

In this section, the MS error resulting from the error between the true wavefront tilt, 0 t , 
and the wavefront tilt “seen” by the subaperture is computed using a method similar to that of 
Yura and Tavis [24]. The Hartmann sensor averages the wavefront tilt, V<j>, over the subaperture. 
In doing so, the phase slope is approximated in a piecewise linear fashion. As can be seen from 
Figure 3.3(a), this is a good approximation when r, is large. However, as r„ decreases, the averaged 
wavefront tilt diverges from the true wavefront tilt, as can be seen from Figure 3.3(b). 
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Figure 3.3. Hartmann Sensor with (a) large r 0 and (b) small r, 


The MS error between the z-tilt and the estimated tilt can be expressed as [24]: 

({St-tf) (3-25) 

where 0 , is the true wavefront tilt, called z-tilt, over the entire aperture and 0 , is the estimate of 
this tilt. This estimate, 9 t , is assumed to be equal to the instantaneous angular position of the 
centroid, 9 e , which is equivalent to the a in Figure 3.2. 0 e is commonly called the 9 -tilt. Expanding 
Equation 3.25, we obtain: 

( (Ot - Oc )^ - 20, • 0 e + 0?) (3.26) 

Recall that 0, and 9 e were defined in Section 2.3. These definitions will be used to evaluate each of 
these three terms individually. 

Using Equations 2.39 , 2.40, 2.44 and 2.45, the cross correlation between 9 -tilt and z-tilt, 
(Ox 0c) can be expressed as: 

(ft -9c) = (*•* + #I ■ y)W{z,y)4(x,y)dzdy ■ J W(z,y)V<l>(x,y)dxdy) (3.27) 

Using integration by parts, this becomes: 

( 0 , - 0 e) = (^Jl (xix + yiy) W{xi,yi)4(xi,y{)dx\dyi 

1 V7W(z2,y 3 )<t>(z2,y2)dx2dy2'j S ) (3.28) 





Using the definition of the gradient operator in Equation A.3, performing the dot product and 
simplifying yields: 


(*» Oc) = 

X [* (*’ ** 0 “ 6 ( Xj + j)] " ct (“f ) d*i<l*J d lhdya^ (3 29) 

since all terms in Equation 3.29 are deterministic except the expected value operator can be 
taken inside the integrals to yield: 

x |$ ^x a - 0 - 6 (x-i + ^j rect diKixjdyidyj (3.30) 

An expression for the correlation function, (^(*i,yi)^(xa,!fe)), is needed to continue. It is shown 
in Appendix A.l, that, using the Kolmogorov structure function as described by Wallner [21:pg. 
1774]: 

it 


D (x, x‘ ) = 6. 


88 


(*-*') 


(3.31) 


the correlation function can be expressed as: 


(<t>(xi,vi)W**y2)) = ■ 1 ^ + ^ ya +r*(o,o) (3.32) 


Substitution of Equation 3.32 into Equation 3.30 yields: 


_i 6 88 ^ V(»a-*i ) 3 + (ya-yi) a j 5/3 

[<5 ^2 - 0 - 6 (x 2 + 0J rect dx!<ix a dyidife (3.33) 


<* ■*> = N 6 - 88 ( — ■■ j 


Using the sifting property of delta functions to integrate over x a , making the variable substitutions 
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given by Equations A.17 and A.18 and using Fourier Transform analysis, Equation 3.33 becomes: 


<«■*>-£ 


6.88 


2rV 3 







(3.34) 


Substituting in for the definitions of the tri and rect functions, setting the limits on the integrals 
using these definitions and performing numerical integration, Equation 3.34 can be expressed as: 


(£•£) = 


13.0934 

i 2 r! /3 d»/3 


(3.35) 


Recalling the definition of 9 t , Equation 2.39, and also recalling Equation 2.44, the mean 
square value of the r-tilt, < 9^ >, can be written as: 


<*>-»)'((/ 


(xx + yy) ^(x, y)W(x, y)dxdy 


)') 


(3.36) 


Using a very similar analysis as that used to determine ^ ■ 9 C ^, it is shown in Appendix A.2 that 
Equation 3.36, becomes: 

13.6183 


w- 


k 2 rV 3 d }! 3 


(3.37) 


Likewise, the MS value of 0 -tilt, < 9% > can be expressed using Equations 2.40 and 2.45 as: 

<« 




•,y)W{x,y)dxdy S j ^ 


(3.38) 


Equation 3.38 is also evaluated in Appendix A.3, using a similar approach, and becomes: 


(«)- 


12.803 

fc 2* /3 <fl/3 


(3.39) 
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3.2.1 Combining < 0$ >, < > and <0t -9t > to Find Eg t , Substituting Equations 3.35, 

3.37 and 3.39 into Equation 3.26, will result in: 



0.235 

fc 2 r* /3 di/3 


(3.40) 


which has the same functional dependence as a similar quantity determined by Yura and Tavis [24] 
for a circular aperture. Equation 3.40 can also be written as: 



0.235 




(3.41) 


Equation 3.41 has units of dimensionless slope (meters over meters). However, the (fcr 0 ) 2 term has 
units of: 


( radians \ 2 / meters \ 2 / radians \ 2 

meter ) \ r o cell / \r 0 cell ) 


(3.42) 


Hence, the product, < (Q t — &c) 2 > (kr 0 ) 2 will have units of (radians/r 0 cell) 2 , which is desired. 
Therefore, define: 


(3,3, 

3.2.2 Total Hartmann Mean Square Error Combining Equations 3.24 and 3.43 will give the 
total MS error for the Hartmann Sensor due to both shot noise and the error between (/-tilt and 
z-tilt. 

Eb = [ ^ V " + <°- 37 > ! W’] + °- 235 (£)'* f » £ s 1 (3M) 

\ 2!2^[^ + (0.37) j ]+0.235 (a)"*(i^^) 2 for * > 1 

Some observations regarding Equation 3.44 are in order. First, this result shows that some error 
is inherent, even for “perfect” seeing conditions, i.e. d/r 0 < 1. That error is determined by the 
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amount of light, N, and the diffraction limited spot size in the detector plane. Second, the shot 
noise portion of the error is strongly dependent on the angular size of the source, 0. Hence, as 
the source becomes larger, i.e. an artificial guide star as opposed to a natural star, the error will 
increase. And third, the error is dependent on the seeing conditions, given by the factor d/r„. 
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IV. Shearing Interferometer Wavefront Sensor 


As for the Hartmann Wavefront Sensor, two major sources of error will be examined for the 
Shearing Interferometer. Those two sources can be described by the error between s-tilt and r-tilt 
in the wavefront measurement and the shot noise effects. These two effects are independent and 
hence can simply be added to give the total mean square (MS) error. For notation simplicity, the 
total MS error for the Shearing Interferometer will be defined as: 

Es = E Sp + Estt (4.1) 

Es P is the MS error due to the photon, or shot, noise effects. Es, t is the MS error due to the 
difference between the true wavefront tilt, and the tilt “seen” by the Shearing Interferometer. A 
similar argument as given in Section 3.2 applies. To match units with the Hartmann analysis, 
Es will be manipulated to have units of (radians/r 0 cell) 2 . An expression for Es p is derived in 
Section 4.1. Likewise, an expression for Es,, is derived in Section 4.2. 

4.1 Shot Noise, Es p 

It was shown in Chapter II that the count vector statistics for our AC lateral Shearing 
Interferometer model match those of Goodman’s DFT Fringe Estimator. Thus, the result derived 
by Goodman [8] for the RMS error in estimating due to shot noise can be used for our model 
as well: 

= 

where a is the quantum efficiency of the detector, A is the area of the detector element, r is 
the integration time of the measurement, is the intensity of each sheared beam and V is the 
visibility. The notation A^|* indicates this is the RMS error for the x shear section of the Shearing 
Interferometer. 
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It is important to relate Ii and In to the incident beam intensity. The light coining into the 
Shearing Interferometer must be split at least twice. The first split divides the light between the 
x and y shear sections of the Shearing Interferometer. The second beam splitter divides the light 
into two beams which are sheared and interfered. Assuming 50-50 beam splitters, the individual 
beam intensities can be described as: 


T r I T . I 

h - h = ~ — I\ + h = j 


(4.3) 


Since the two beam intensities are equal, the visibility is equal to the complex coherence factor [8:pg. 
163]. The complex coherence factor for a Gaussian source of RMS width <r, was found in Chapter III 
to be: 


/ii 2 (Ai, Ay) = exp 



(Ax 2 + Ay 2 ) 


(4.4) 


where Ax and Ay, are the separations of the points in the x and y direction. Therefore Equation 4.4 
can be used to write the visibility as a function of d, the aperture size, and a = Ai/r„, the 
normalized shear. 


V(s,r„) = exp 



(4.5) 


where A is the mean optical wavelength and z is the distance from the source. Substituting Equa¬ 
tions 4.3 and 4.5 into Equation 4.2 yields: 




-t 


aArNI 


exp [-2(5^) 2 (sr„) 2 


(4.6) 


which has units of radians. To effectively compare this with the results for the Hartmann sensor, 
the units on Equation 4.6 need to be converted to (radians/r 0 cell) 2 . First, Equation 4.6 needs to be 
converted to a slope. This is done by dividing Equation 4.6 by Ax, which gives the phase change, 
(in radians) per amount of shear, (in meters). Further, squaring Equation 4.6, to make it a MS 
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error instead of a RMS error, produces: 


*!»!,(*' r °) 

Ax 2 




where N to t = aArNI, which represents the total average photon count over the aperture for each 
integration time. N to t is comparable to N used for the Hartmann sensor. Equation 4.7 has units 
of (radians/meter) 2 . Thus, by dividing by 1/r 2 (or effectively multiplying by r 2 ), which has units 
of (meters/r 0 cell) 2 , the desired units will be obtained: 


x(*- r ») r o 


4r 2 

Ntot&x 2 


-[<(£)'<->■] (a?)’ 


In Equations 4.6-4.8, the |x designates that this is the error in the x direction. If we assume that an 
identical error will occur in the y shear section of the Shearing Interferometer, then the total er or 
will be twice that indicated in Equation 4.8. Further, recalling that the s = Ax/r„ Equation 4.8 


becomes: 


r «) n 4 


= 2 jw exp 


[«(S)W 


Thus, the MS error due to the shot noise can now be defined as: 


Es P = 


_°i*(*> r o) _ 8 

S 2 ~ Nt 0 ts 2 


- [« (£)’<->’] (^)’ <«■> 


The terms in the exponent can be rearranged to yield: 


Es ' = 


(4.11) 


Recalling Equation 3.23, Equation 4.11 becomes: 


Es p = 


8 

N tot s 2 




(4.12) 
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4-2 Mean Square Error between ;-tilt and s-ttlt, Es,, 

This section examines the second of the two major sources of error in the Shearing Interfer¬ 
ometer. As for the Hartmann Wavefront Sensor, we are interested in the MS error between the 
true tilt and the tilt estimate given by the sensor: 

(4-13) 

where 9 t is the true wavefront tilt over the aperture and 9 t is its estimate. In the Hartmann Sensor, 
it was assumed that 8 t = 9 C , where 9 C is the instantaneous angular position of the centroid in the 
Hartmann subaperture. A similar approach will be taken for the Shearing Interferometer. The 
estimate of the wavefront tilt, 9 t will be assumed to be equivalent to 8, which is the tilt measured 
by the Shearing Interferometer. Expanding Equation 4.13 then yields: 

( (*‘ “ *’Y) = (# - 2 *‘ ’ + 3) (4T4) 

The mean square value of was evaluated in Section 3.2 and found to be: 



13.6183 

k 2 rV Z d l l 3 


(4.15) 


The terms in Equation 4.14 that remain to be evaluated are the mean square value of s-tilt, < 8% >, 
and the cross-correlation between z-tilt and s-tilt, < 9 t ■ 9, >. 

Using Equations 2.41 and 2.46, the mean square value of s-tilt, < 9% >, can be expressed as: 
( P ') = ({h) (*/ A x ^ W{x, y)didy A y Jf ~ W(x, y)dxdy) ^ (4.16) 
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It is shown in Appendix B.2 that evaluating Equation 4.16 yields: 


(fi) = Jiff (A*(*i,yi,Ax)A*(x 2 ,9a,Ax)) W(xi,yi)^(*s,yi)d*idyida: 2 dtfj (4.17) 

An expression for (A^(xi,yi, Ax)A^(x 2l y 2l Ax)) is therefore needed. Recalling the definition for 
A <(>: 

A^>(x, y, Ax) = ^(* “ y) ~ <^(x + y) (4.18) 

and using the structure function given by Equation 3.32, Equation 4.18 becomes: 

(*?) = k 2 d 4^ 2 //3 if [>/ ( * 2 ~ Xl + A * )2 + ^ ~ **) 2 + >/(**“*i-A*) a + (»j-yi) 2 

_ 

-2y(x 2 —xi) 2 + (y 2 - yi) 2 

xrect rect rect ^ j rect dx 1 dy 1 dx 2 dy 2 (4.19) 


A variable substitution can be made so that the integral in Equation 4.19 will be in dimensionless 
form. This includes a substitution of s = Ax/r„, which is defined as the normalized shear. Using 
these substitutions and Fourier Transform analysis, Appendix B shows that Equation 4.19 becomes: 


a\ 0 13544A(*fr) 

**d»/»r5 /3 


where /„ (s 2 ^-) is defined to be: 


fv (X) = (x)- 2 J* J\ 2 - |A*|) (2 - |A^|) |((A* + 2x) 2 + A^> 2 | * 
+ |((A4 - 2x) 2 + A^ 2 | f - 2 |(A< 2 + AV- 2 ! 1 


dA£dAV> 


(4.21) 


and is plotted in Figure 4.1. Equation 4.20 is the simplest analytical solution than can be obtained. 
Note that any numerical solutions of Equation 4.21 must be for a given value of jp~- Also, note 
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that from Equation 2.47 



(4.22) 


Hence, Equation 4.20 will be solved for a family of values from, 0 < < 1/2, as it is 

unreasonable to assume that the wavefronts would be sheared more than one half the subaperture 
diameter. 

Attention is now turned to the middle term of Equation 4.14, (o t By use of Equa¬ 

tions 2.39, 2.41, 2.44 and 2.46, this can be expressed as: 




f x<fi(x,y)W(x,y) dxdy 

kd* 

12 


+ y 


f y<t>(x,v)W(x,y) dxdy 
kd* 

12 


/ A0(z, y, Ax)W(x, y)dxdy 
Axkd 2 


+ ii 


/A^(z,y,Ay)iy(x,y)dgdy j \ 

Ay/M 3 J/ ' ’ 


Using a very similar analysis as that used to determine (fy, Equation 4.23 is evaluated in Ap¬ 
pendix B.3 and is found to be: 



0.406325/c (sfr) 

*2 rf l/3 r 5/3 


(4.24) 
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32.2V 



where f e (s^) is defined as: 


U (*) = *“ l f f Mil- \M\) (2 - |A*|) 

X |j(A£ + x) 2 + Atf 2 |*-|(A£-z) 2 + A<l> 2 


dA£dAV> 


(4.25) 


and is plotted in Figure 4.2. Equation 4.24 is the simplest analytical form of an expression for the 
cross correlation between 3-tilt and s-tilt. To proceed from here, a value for jp~, the normalized 
shear, must be specified and Equation 4.25 solved numerically. 


4.2.1 Mean Square Error Between, z-iilt and s-tilt, Es„ Recalling Equation 4.14, the total 
MS error due to the difference between 2 -tilt and s-tilt can be obtained by combining Equa¬ 
tions 4.15, 4.20 and 4.24. Doing so yields: 



13.6183 _ 2 • 0.406325/e (sfr) 0.13544/, (sfr) 

F 2 r* /3 dl/3 i 2 rfl/3 r 5/3 + *2rfl/3 r 5/3 


(4.26) 
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where /„ (s^f) and /<.(s^) are defined by Equations 4.21 and 4.25 respectively. Note that Equa¬ 


tion 4.26 can be re-wiitten as: 

^(<f« ~£) 2 ) = [13.6183-0.8126/ e (s^) +0.1354/, (s^)j ^ (4.27) 

which gives the MS error in terms of slope, (meters/metere). Therefore, as in Section 3.2.1, to 
convert the units to (radians/r, cell) 3 , a multiplication of ( kr 0 ) 3 is required. Thus: 

Es., 3 ((», = [13.6183-0.8126/. (.i) + 0.1354/. (.$)] (£) 4 (-£==)’ 

(4.28) 


4.3 Total Shearing Mean Square Error 

Using Equations 4.1, 4.12 and 4.28, the total MS error for the Shearing Interferometer can 
be described as: 

= a!CP + [13 6183 - 0.8126/. ( 5 ^) + 0.1354/. (.^)] (i) ‘ (^f)’ 

(4.29) 

Some comments about this result are in order. First, as in Equation 3.44, Equation 4.29 shows 
that some error is inherent even for “perfect” seeing conditions, d/r 0 <C 1. This minimum error is 
determined by the total photon count, Ntot , and the shear, s. Figures 4.1 and 4.2 show that the 
f c (s r ~f ) and f v {s^f) change very little over the range of values we are interested in. Thus 
the MS error between s-tilt and z-tilt is only weakly dependent on the shear. 
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V. Comparison of Hartmann and Shearing Wavefront Sensors 


The total error variance for the Hartmann Wavefront Sensor, found in Chapter III, is: 


( 2^ [/J> + (0.37) 1 (^)’l +0 235 (A)'‘ )’ for ± < 1 

r> rr —_ f 

| 2^ [/9 J + (0.37) 2 ] + 0.235 (£) “ * (l^f ) * for £ > 1 

Likewise, the total error variance for the Shearing Interferometer, found in Chapter IV, is: 


E s = K 2 *)^* 2 ] + [13-6183-0.8126/ e (s^) + 0.1354/, (sy)] (yj 


~ ^ ( radians \ 2 
cefiy 




(5.2) 


Recall that the N in Equations 5.1 and 5.2 represents the total number of photon events within the 
subaperture during the integration time. The ratio of the angular size of the source to the angular 
resolution of the optical system, limited by the atmosphere, was defined in Equation 3.23 to be: 



(5.3) 


The normalized shear, s, was defined as the ratio of the absolute shear, Az, to the atmospheric 
correlation length, r„. 


s = 


Az 

r 0 


(5-4) 


The functions f„ (s^) and f e (s^f) are given by Equations 4.21 and 4.25 respectively and are 
plotted in Figures 4.1 and 4.2, respectively. And, finally, d/r 0 is the ratio of the subaperture 
diameter to the atmospheric correlation length. 


Several items about Equations 5.1 and 5.2 should be noted. First, some error is inherent in 
wavefront measurements, even for “perfect” seeing conditions, i.e. d/r 0 < 1. This minimum error is 
determined by the total photon count. N, and the source size, j3. Second, examination of Figure 5.8 
reveals that the term [13.6183 - 0.8126/. (s^f) + 0.1354/, (s^)] is approximately equal to 0.235 
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for most values of of interest. Hence, the error due to the difference between j-z-tilt and s-z- 
tilt, given by the second term in each equation, is relatively the same for both wavefront sensors. 
Third, due to the source sue term in the exponent of Equation 5.2, the Shearing Interferometer is 
much more sensitive to changes in the source size than is the Hartmann Wavefront sensor. 

The total RMS error for the Hartmann Wavefront Sensor and Shearing Interferometer is 
plotted in Figure 5.1 versus d/r 0 for -V = 2; s = 0.1, 0.2, 0.3, 0.4 and 0.5; and 0 = 0.01, 0.10, 0.50 
and 1.00. Looking first at the 0 — 0.01 plot, recall that this source size represents a very small 
(almost a point) source. Notice that the Hartmann Sensor performance improves significantly as 
d/r 0 increases, for d/r 0 < 1 and then flattens out for further increases in d/r 0 . This results from 
the spot size iu the detector plane of the Hartmann being inversely proportional to the subaperture 
diameter. As the spot size decreases, the Hartmann sensor determines the centroid more accurately, 
and hence the shot noise decreases. However, for d/r 0 > 1, the spot size is limited by r 0 instead of d 
and no further reduction in the spot size occurs. We also notice that the RMS error for the Shearing 
Interferometer decreases with increasing shear. As the shear increases, the phase difference will also 
increase while the noise in the measurement will remain relatively the same. Hence, the signal to 
noise ratio will increase causing the shot noise portion of the error to decrease. This is represented 
by the s~ 2 in the first term of Equation 5.2. 

Now, looking at the 0 = 0.10 plot, we see the same trends observed for the 0 = 0.01 plot. 
However, close examination shows that all the curves have been shifted up slightly, indicating 
that the RMS error increases for increasing source size. This is more easily seen in the plots for 
0 = 0.5 and 1.00. However, in these last two plots, we also notice that the consistent decrease in 
RMS error for increasing shear no longer occurs. Instead, error decreases until some “optimum’’ 
shear is reached and then begins to increase again. This results from the Shearing Interferometer’s 
sensitivity to the fringe visibility, which is inversely proportional to the size of the source. 

The total RMS error is again plotted in Figures 5.2 and 5.3 for the same conditions as above 
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except for N = 5 and 10, respectively. The same trends as mentioned in the above paragraphs are 
noticed. However, we see that the total RMS error decreases with increasing photon count since 
the shot noise error is inversely proportional to N. 


Figure 5.4 is a plot of the RMS error for N — oo. This plot represents the portion of the 
total RMS error due only to the difference between g-tilt and z-tilt for the Hartmann Sensor and 
the difference between s-tilt and z-tilt for the Shearing Interferometer. We can clearly see that, in 
the absence of shot noise effects, the two sensors show nearly identical performance. 

Likewise, Figures 5.5-5." represent only the shot noise portion of the total error for the 
Hartmann and Shearing Interferometer-given by the first term in Equations 5.1 and 5.2 respectively. 
Comparing Figure 5.5 with Figure 5.1 we see that they are nearly identical, indicating that shot 
noise effects are the dominate noise source for both the Hartmann Wavefront Sensor and the 
Shearing Interferometer, for the photon count values examined. This same conclusion can be 
drawn by comparing Figures 5.6 and 5.2; and Figures 5.7 and 5.3. 


As mentioned above, there seems to be some “optimum” shear, at which the RMS error for 
the Shearing Interferometer is a minimum. Since Figures 5.5-5.7 indicate that the majority of the 
error results from the shot noise component, that term will be examined to determine the optimum 
shear. Recalling Equation 4.12, the shot noise component of the total Shearing MS error; taking 
its derivative with respect to s; setting it equal to zero and solving for s gives the optimum shear 
as a function of source size: 


1 

Sop, ~2;70 


(5.5) 


Hence, as expected, the optimum shear is inversely proportional to the source size. Equation 5.5 
is plotted in Figure 5.9. 
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Figure 5.1. RMS Error vs d/r a for /V = 2 
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Figure 5.2. RMS Error vs d/r 0 for N = 5 
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Figure 5.4. RMS Error vs d/r 0 for N 
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Figure 5.9. Plot of Optimum Shear, s, vs. Source Size, /?. 
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VI. Conclusions and Recommendations 


6.1 Conclusions 

This research has conducted an analysis of the mean square error resulting from wavefront 
measurements obtained from the Hartmann Wavefront Sensor and the Shearing Interferometer 
when viewing extended objects through a turbulent atmosphere. Two sources of error were con¬ 
sidered: the shot noise error and the j-z-tilt and s-z-tiit error. The shot noise error results from 
inaccuracies in determining the spot centroid, for the Hartmann Sensor, and the phase of the fringe 
pattern, for the Shearing Interferometer. Likewise, the y-z-tilt error and the s-z-tilt error results 
from the difference between the true wavefront tilt and the averaged wavefront tilt over the finite 
sized subaperture. 

Results indicated that the major contributor to the total error, for the photon counts ex¬ 
amined, was the shot noise effects. The shot noise error had a strong depend-nee on the source 
size, in addition to dependencies on the photon count and amount of shear, (for the Shearing In¬ 
terferometer). When the shot noise error was not considered, the two detectors performed nearly 
identically. 

However, for the more realistic-photon limited-cases, the Hartmann Wavefront sensor showed 
much better error performance than the Shearing Interferometer. The only exceptions were for small 
source sizes and small d/r 0 . 

6.2 Recommendations for Further Study 

The model used for the Shearing Interferometer encoded the wavefront tilt in the phase of 
an ac modulated fringe pattern. Goodman’s DFT Fringe Estimator was then used to extract the 
phase from this fringe pattern. The results indicate that the major error resulted from this portion 
of the wavefront detection process. However, there are versions of the Shearing Interferometer 
that do not use ac modulation and fringe detection [11], It is therefore recommended that the 
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analysis presented in Chapter IV be repeated for one of these other configurations, without the 
ac modulation. This would establish whether the unusually large error is a result of the shearing 
process or the ac modulation and fringe detection process. 
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Appendix A. Derivation Details for Hartmann Wavefront Sensor 


This appendix provides details of the derivations presented in Chapter III. 

A.l Derivation of Cross Correlation between g-titt and z-Tilt, (§ t • 

From Equation 3.27 

(h • = pp l^jj (xx + y-y)W(z, y)4>(x, y)dxdy • JJ W(z, y)V<f>{x, y)dardy^ (A.l) 

This second integral is very similar to one evaluated by Welsh and Gardner [23:pg. 1914] using 
integration by parts. This approach yields: 

(ft £) = pp (^jj {xix + yiy)W(x u yi)<Hx u yi)dxidyi 

I V2 ^* 2,w M* 2, ^ dl2dy2 )) ( A - 2 ) 

Note that in Equation A.2, dummy variables, denoted by subscripts 1 or 2, are used to indicate 
the separate integrals. Further, note that V 2 indicates that the gradient operator only operates on 
functions of x 2 and 1/2 

Recalling the definition of the subaperture, W{x,y) (Equation 2.42), and recalling the defi¬ 
nition of the gradient operator: 

v=i h +i h (A - 3 > 

T 7W(z,y) becomes: 

= h h‘ (I) " ct ©l i (3)'“*(?)]» 

= [-<(*-|) +« (* + 0] reel (J) i 

+ ("-1) +*(»+|)] r “ , (|)w ( A “) 
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Substituting this relationship into Equation A.2 yields: 


■ (- jf <K*7, va) {[s (*a + 5 ) -' 6 (««- 0 ] «ct ( 5 ) i 

+ [ 6 (» + 0 " 6 O'* ~ O] reCt (?) »} dijd »>)) 


(A.5) 


Separating each integral in Equation A.5 into its z and y components and performing the dot 
product results in: 


(ft e c ) = p^^^2:iW(zi,yi)^(*i,yx)dxidyi 

x jf <£(z 2 ,y 3 ) ^z 2 - 0 (* 2 + j)] rect (y) dzjdyj 

+ ff yiW(zi,yi)^(x u yi)dxidyi 

* jf <K*a. Vi) |$ (ife - 0 - 4 (»2 + 0 j «*t (y) dz 2 dj0 (A.6) 


Note that all actions of the first two integral pairs will be duplicated on the second two integral 
pairs. Thus, the second two integral pairs are dropped and the entire expression is multiplied by 
two to account for them. This yields: 


(ft -9^ = *1 W(*i> ,yi)d*idyr 

* Jl <H x i, W) Jf (*a - |) ~ 6 (*j + 0] rect dx 2 dy0 (A.7) 
Combining the integrals in Equation A.7, produces: 

(Off>c) = ^(JUI x i w ( x i>yi)H x i>yi)<l>( x 2iy2) 

*’-£H 



^z 2 + 0 rect {fjj dzidz 2 dyidy0 (A.8) 
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Note that all terms in Equation A.S are deterministic except Thus, the expected value operator 
can be taken inside the integrals. This results in: 


(f* £) = pp J[jf *i^(*x.yi)^(*i.yi)^{*j,yj)> 

X ^ ^x 2 - 0 “ 6 + !)j wet (^) d*idxjdyidyj (A.9) 

An expression for the correlation function, {^(*i,yi)^(*j, jft)), is needed to continue. However, 
in atmospheric optics, it is more common to use structure functions rather than correlation func¬ 
tions [22]. The correlation functions and structure functions are related by: 




(A.10) 


where r*(0,0) = j (^ 2 (*i, yi ))+5 (o 2 (* 2 , !&)) and is constant since the phase is assumed to be wide 
sense stationary [8:pg. 376]. The structure function, D+, has been determined by Wallner [21:pg. 


1774] to be: 




l ^ ( £ ~^) 

D ^x, x'J = 6.88 ' 


(All) 


Therefore, the correlation function can be expressed as: 


(^(*i>yiM*2,j/2)) = -- 


6.88 ( v(* 2 -*i) 


2 + (ift-yi) 2 ^ 


+ r # (o,o) 


(A-12) 


Substituting Equation A.12 into Equation A.9, produces: 


x ^ 6 ^x 2 - 0 - 6 (*2 + j wet (^) dxidxjdyidyj (A.13) 


Note, if f* is separated out from the rest of the integrals and the limits imposed by the subaperture 
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are set, the resulting integral will be a constant multiplied by an odd function integrated from 
— j to j and hence will become zero. Therefore, the T* term can be dropped. Recalling the 
definition of the subaperture function, Equation 2.42, Equation A.13 becomes: 


(*'' *) = *37 (7) ( d ) [' (* i - 5) -■>(*=+ 1 )] (7) 


6.88 ( 

2 V 


y/(x 2 -jr t ) 2 4-(P3-yi) 2 j 


5/3' 


dxjdxjdyidy? (A.14) 


Using the sifting property of delta functions to integrate over xj, [6:page 53], 


J f(a)S(a - x 0 )da = f(x a ) 


(A.15) 


Equation A.14 becomes: 


('* •»«) = £*•’“* (7) r “‘ (7)"" (7) 


6.88 ( 1) 2 + {y2-y i) 2 

2 


5/3 


6.88 (>/(—f -*i) 2 + (y2-yx)* 


5/31 


dxidyidyj 


(A.16) 


To proceed, the following variable substitutions are defined: 


Ay = yj - yi 


and 


£y — 


m +yi 


(A.17) 


Thus, 


yj = -Ay+Ey 


and 


yt = Ey - -Ay 


(A.18) 


This choice allows the Jacobian to be one. Equation A.16 then becomes: 
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{’• ’') =J " 11 '" 1 (7) " ct ( Sg j' Ag ) rect (^f^) 


- — - 5/3 . _ 

6.88 I -M- * 1 ? + ( A y)’ ^ 6.88 / */(-$- ti) 1 +(Av) 7 } 

" 2 l " / + 2 l - j 


5 / 3 - 


xdxidAydEy 


Note, only two terms in Equation A. 19 depend on Ey. Thus, for a moment, just look at: 


(A-19) 


/ re "( ^-" ) rect ( Sg / Ag ) 


(A.20) 


Define: 

^Ay + Ey 

Or =: -- 

d 


Then Equation A.20 becomes: 



(A.21) 


d 


J rect(a)rect ^ar - da 


(A.22) 


which, since Ay/d is a constant for this integral, is immediately recognized as the correlation of 
two rect functions. Thus, Equation A.20 can be evaluated by referencing Gaskill [6:page 172] and 
using Fourier transform analysis. 

J rect rect ( ? y - j dEy = d tri j (A.23) 

Substituting the result of Equation A.23 into Equation A.19 results in: 
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<« *> - £/— (?) [' - (?)] [-t 2 


6.88 ( \/(i + x i) i + (^y) 2 


dzidAy 


(A.24) 


This double integral can be broken into two double integrals yielding: 




2 + *i J + (Ay) 2 l dzidAy 


Recalling that: 




2 ~ *i ) + (Ay) 2 l dzidAy (A.25) 


(?)- 


-1 < «i < f 

otherwise 


(A.26) 


m= 


- d < Ay < d 


otherwise 


the limits on the integrals in Equation A.25 can be set as: 


(A.27) 


d d 

- < Zi < — 
2 ~ ~ 2 


— d < Ay < d 


(A.28) 


and Equation A.25 becomes: 




'LL’t*'(*-£*){{*-*') +(i!,)! ) d * ,dA » 


(A.29) 
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Numerical techniques must be used to proceed from this point. However, Equation A.29 will first 
be written in terms of dimensionless quantities to facilitate the numerical integration. To do so, 
let: 

~ <A - 30) 

with limits: 

- 1 < £ < 1 (A.31) 

And let: 

A y / = ^ - dAy / = idA y (A.32) 

a a 

with limits: 

- 1 < Aj/ < 1 (A.33) 


Substituting Equations A.30, A.31, A.32 and A.33 into Equation A.29, produces: 




y/s 


*)«■«' 


(A.34) 


Performing algebraic simplification yields: 


(* »«) = [/_ j_A (1 - la.K'l) (d + *1)’ + (2A # ')')‘"d«idA»' 


- / f (1 - Ay'll ((1 - «;)’ + (2Ay') i ) V ‘dr',dAy' 


(A.35) 


Using numerical integration, the difference of these two integrals is found to be 2.014. Thus the 
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cross-correlation for a square subaperture can then be expressed as: 



13.0934 

fc 2 r * /3 di/3 


(A.36) 


A.2 Derivation of Mean Square Value of z-Titt, < 9% > 

From Equation 3.36: 

(<%) = (^r) ^ i^jj (*■* + *•*) *(*> y)W(z, y)dxdyj \ (A.37) 

Which can be written as: 


• Jf (*2* + y2y)^(*2,y2)W'(* 2 ,y 2 )dz 2 dy 2 ^ (A.38) 

Next, the dot product will be performed which will eliminate the unit vectors. Further, just as 
in Section A.l, the dot product will form two integrals which will be identical except for different 
variables. This duplicate integral can be accounted for by multiplying one of the resulting integrals 
by two and dropping the second. Doing this results in: 

(ff) = -fc ij g- (^Jfff *i*20(*i,yi)^(*2,Jfe) W'(xi,yi)W(*j,y J )dxidyidx 2 dpj^ (A.39) 

Since <t> is the only random process in Equation A.39, the expected value operators can be moved 
inside the integral yielding: 

) = \r]t Jfff x ^ x t{ < t > ( x ^yt) ( t>{xz,y 2 ))W(xi,y x )W{x 2 ,]n)ddxidyidz 2 dy 2 (A.40) 
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In Section A.l, Equation A. 12 , the correlation function, {<f>(zi,yi)<p (12.S/2)). was found to be: 


= tl)a fB ' f(ya — j +r # (o,o) (A.4i) 

Note that if the r^(0,0) term is separated out from tho rest of the integrals and the limits imposed 
by the subaperture are set, the resulting integral will be a constant multiplied by an odd function 
integrated over symmetric limits and hence will become zero. Thus the 1^(0,0) term can dropped 
here as well. Including these substitutions and taking out the constants produces: 

(ff ) = Jiff ((*2 — *i) 2 + (V 2 ~ yi) 2 )*^ 6 W (*i,yi) W (*a>yj)d*idyidx 2 dy 2 

(A.42) 

The subaperture function, W(x, y), serves only to set the limits on the integration. Setting those 
limits yields: 


/ a \ -6.88-12 2 / d/2 f 4/2 [ 4,2 [ d12 „ , 2 ,2x5/6, , . , 

(ff) = ~ r; 5/3 h. / / / / *i* 2 ((* 3 -*i) 2 + U/2-Sh) ) dndyidijdjft 

x ' k-rq'^d* J-dpJ-dpJ-dpJ-dp 

(A.43) 

Equation A.43 can be written in terms of dimensionless quantities to facilitate the numerical inte¬ 
gration. Let: 


c _ *1 

^ d/2 

-* d^i = ^dxj 

- 1 < 6 < 1 


c _ 12 

^~dj2 

— = jdl 2 

- 1 < £2 < 1 


, _ » 
d/2 

— dV>i = ^dyi 

- 1 < tfi < 1 


, _ Vi 
d/2 

— d^ 2 = ^dj /2 

- 1 < l/>2 < 1 

(A.44) 
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Using these variable substitutions and performing some simplification yields: 


)= 262 vw r £ d w L L L L 66 ( i6 ■ 6)2 + { * 2 ~ v ’ i)2 ) ; 

(A.45) 


Using numerical integration to solve this 4-D integral results in: 



13.6183 

|.2 r 5/ 3 ( fl/3 


(A.46) 


A.3 Derivation of Mean Square Value of g-Tilt, < 0% > 
From Equation 3.38: 



V</>(x,y)fV(x,y)dxdy 


(A.47) 


Which becomes: 

(^) = ^ Jf ^i<i>(xi,yi)W{x 1 ,y l )dxidyi J 


V2<H x 2>y2)W(x 2 ,y2)ddx 2 dy2 


) 


(A.48) 


Welsh and Gardner [23:page 1914], have shown that: 


Jl W(x, y)V<t>(x, y)dxdy = - VW(z, y)<i>(x, y)dxdy 


(A.49) 


Therefore, Equation A.48 becomes: 

( 0 ?) = v i^( a; i.yi)^i.yi)diidyi • Jl V 2 W(i2, 2/2)<£(*2i 2/2)d^2dy2^ (A.50) 
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Combining the integrals, results in: 


(^) = (^JIH ^^(xuyi) ■V2W(x2,y7)4>{xi,yi)4>(x2,yi)dxidyidx i dyi^ (A.51) 

Again, since <j> is the only random variable, the expected value operators can be moved inside to 
yield: 

(^) = kh* J[[f ' 7lW ( Xl ' y i)' v 2W , (x 2 ,y2)(^(ii 1 yi)^(i2,lft!))da:idy l dx2d»j (A.52) 

Equation A. 12 transforms the correlation function to a structure function, resulting in: 

(*1)= Jjff v 1 ^(* I , w )-v a wf(* a ,»»)(r # (o,o) 

(A . M) 

Recalling Equation A.4, which expands out VW(x,y), Equation A.53 can be written as: 

x {(K I1 + 0‘ < ( I,_ O]" c ‘(i') i+ K , ‘ + l)- 4 (" - i)]" ct (T) s ) 

■ ([*(* + 1) - 4 (*>- 5)] ■** (?) *-+ K” + 1) - ‘ (»-!)] “ t (t) *) } 

xdx x dyidx2dy2 (A.54) 

Performing the dot product yields: 
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«{[«(*.+0 -< (*. - 1 )] [■* (*>+ i ) -• (*> - 1 )] ■"* ( 7 )( 7 ) 

+ [«(»>+ 5 )-*(«-!)] [<(» + 0 -»(«- 0 ]»ct(y)-ect( 2 )} 
xdxidyxdx2<iy2 (A.55) 

Examination of Equation A.55 reveals that it can be split into two 4-D integrals differing only by 
what is an 1 in one delta or rect function is a y in the other. Thus, as before, the second integral 
can be dropped and accounted for by multiply the remaining integral by two, yielding: 



x { K*‘ + 0 " ‘ ” 0] [‘ (*’ + 0 - 5 (* ! - 0] recl (?)"“ (?)} 

xdx!dr/idx 2 dj /2 (A.56) 


Multiplying out the delta functions will result in: 



To make the next few steps easier to visualize, Equation A.57 can be written as: 
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<«> ■ ( 


k 2 d 4 



«ct (|) «ct (£) r*(0,0) - S,S 

* s (* l + f) ^ ( I! + 0 diidyidudya 

- f - (5) - (I) [™ »> - t 5 ( ^EinsEH^ j 5/3 

*5 + 0 6 ^x 2 - 0 dxidyidxjdjft 

~M ,cct ( t )'"'* ( 7 ) 

xd 6 ^r 2 + 0 d*idyidi 3 dyj 

+f - (I) - & [™« - ^ 

xt -05 - 0 dxidyxdx 2 dy 2 | 


r #( o, o) - 6 -f ^ ^ZjjglSz ag^ 


(A.58) 


It can be observed that if the integrals were further split to separate the T functions, all integrals 
with r* functions would subtract out. Hence, we can drop those terms. Now using the sifting 
property of delta functions to integrate over xi and x 2 and recalling Equation A.15, we get: 

^ = (i^){/" ct ® ,ea (7)((-5'H)) 

-/ rec *(7) rect (?) ((!■(■!)) + <»-w) 1 ) 

'Jf rect (7) rect (?) +(w_yi)2 ) dyidyj 

+ Jf rect (y) rect (y) ^(f “ 0 + (v* ~ d yi d V 2 1 (A.59) 

Note that in integral pairs I and 4, the | terms cancel out, where as in integral pairs 2 and 3, the 
| terms add to form d. Thus integral pairs 1 and 4; and 2 and 3 are identical. Once again, the 
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duplicate integral pairs will be deleted and tbe remaining integrals will be multiplied by two to 
account for them. Further, the rect functions serve only to set limits on the integration: 




5/6 


dyidya (A.60) 


Equation A.60 can be written in terms of dimensionless quantities with the following variable 
changes: 


4>i = 


y i 
d/2 


— dyi = -d^>i 


’/'I = ~jj2 ~ d W2 = 2 d ^2 


-1 < V >1 < 1 

- 1 < v>2 < 1 


(A.61) 


which after simplification yields: 


<«>- 


6.88 


2 • 2 5 /3Jb2 r ’/ 3 ^/3 


) [J-i f i ( 4 + ^ ~ V '° 2 ^ 5/6 " ** ~ ^ l|5/3 d * l<W *] (A ' 62) 


Using numerical integration to solve this 2-D integral results in: 


/p\ = J*™ 

N ' k 3 rl ,3 dV 3 


(A.63) 
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Appendix B. Derivation Details for Shearing Interferometer 


This appendix provides details of the derivations for the Shearing Interferometer. 


B.l Derivation of Count Vector, K(n) 


From Equation 2.29: 


K{n) = 2 la Jjj j^l + Vcos + A<p(x, y, Ax)j j IV^x.yJWt^dtdxdy (B.l) 


which can be separated to obtain: 


K(n )= 2 la 




y)W,(t)dtdxdy 


' JH cos + A<l>(x,y,Ax)j W,(x,y)Wt(t)dtdxdy 


The first triple integral is trivial and simplifies to: 


K(n) = 2 la Ar+vJH cos + A*(x, y, Ax)) W,(x, y)IV,(t)dtdxdy] (B.3) 


where A = d 2 has been defined as the area of the detector element. Now, using the relationships: 


t> a + e~ ,a 


e (a+») _ e a e » 


Equation B.3 becomes: 


K{n) = 2/a At + ~ jjj + w,(x ,y)W t (t)dtdxdy (B.5) 
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which becomes: 


A'(n) = 2I(\at + y j e»Vwi(l)d< fl e>*«‘*^W,(x,y)dzdy 

+ | J e-**Pw t (t)dt JJ e-' A * s »^ ) W,(z,y)d*dy (B.6) 

For a moment, just look at the term: 

Jl e i±«*.>.&*)w,{x,y)dzdy (B.7) 

which represents the component of the measurement proportional to the change in wavefront phase. 
Assuming that Ax, the shear, will be adjusted such that A^(x, y, Ax) <C 2x, which will be required 
for accuracy, then the approximation, 

e 7 * « 1 + jx for small s (B. 8 ) 

can be used and Equation B.7 becomes: 

jj W,(x,y)dxdy + j JJ A^(x,y, Ax)W,(x,y)dxdy (B.9) 

Defining: 

AA0(Ax) = Jl A^(x, y, &x)W,(x, y)dxdy (B.10) 

and using the definition in Equation B.10, Equation B.9 becomes: 

A — jAA<l>(Ax) = A (l + ;A^(Az)) (B. 11) 
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And using Equation B.8 in reverse. Equation B-7 becomes: 


eJ*«.**'**)W,(x,y)dz dy = Ae> A «^ 


(B. 12) 


Substituting Equation B.12 into Equation B.6 the count vector becomes: 


7f(n) = 2 la [at + A^ J V W x {t)dt + A^ j e-' A « A *>e-'^W',(t)dt| (B.13) 


which can be recombined to form: 


tf(n) = 2IaA |r + V J cos + A*(Ax) j W,(/)dt] (B.14) 


Recalling Equation 2.28, the limits in the integral can be set as: 


K(n) = 21 a A jr + V jf * * cos (y- + A*(A*)^ dt 


(B15) 


Performing the integration results in: 


K{n) = 2IaA r + ^ sin - + y + A^(Ax)^ - sin - y + A<*(Ax)^ (B.16) 


Using the relationship: 


sin a — sin 0 = 2 cos ^(a + /?) sin ^(a — 0) 

It It 


Equation B.16 simplifies to: 


(B.17) 


tf(n) = 2 IaA jr + ^ sin (y) cos (yp + A<^(Ax)jj (B.18) 






Assuming that the sampling rate is such that: 


XT 

T <1 


^ T 

T < — 

T 


(B.19) 


(note that Nyquist criteria requires r < T/2) then the small angle approximation can be used for 


sin(x) « x for x < 1 


(B.20) 


and Equation B. 18 simplifies further to: 


K(n) = 2IaA r + Vrcos + A^(Ax)^j 


(B.21) 


By using the relationship: 


tN = Po T 


(B.22) 


where r is the integration time, N is the number of samples in the count vector, K(n), i.e. 0 < n < 
N - 1, p 0 is the number of fringes embraced by the count vector and T is the period of the fringe. 
Then Equation B.21 can be written as: 


K( n) = 2IarA 1 + V cos 


( 2xp 0 n 

~ir 


+ 2Ad(Ax) 


(B.23) 


B.2 Mean Square Value of < 0% > 


Using Equations 2.41 and 2.46, the mean square value of s-tilt can be expressed as: 
fi) = ({h) (*/ ~ Ax ,AX) w(a; ’ y)dgdy + ^/ - {X £ y " X) W(*,y)d*dv) ^ (B.24) 


B-4 









Which simplifies to yield: 


W-(ro jl A^(x,y,Ax)W(x,y)dxdy • A*(x,y, Ax)W(x,y)dxdy 

+ kd ?A y Jf A ^( a; -J'« A J') Wr ( ;r -l') dxd y Jf A^(x,y,Ay)W(x,y)dxdy^(B.25) 

Examination of Equation B.2S reveals that the two added terms are essentially the same. Therefore, 
the second term can be dropped and accounted for by multiplying the first by two. Doing this and 
combining the four resulting integrals yields: 

(ff) = ( fc 2^£ ~ 2 f[ff A^(*i.yi, Ax)A0(x2, ya, Ax)W(xi,yi)W(xj,y 2 )dxidyidx 2 dy 2 ^ (B.26) 

Since W(x,y) is deterministic, A<t> is the only random term. Hence, the expected value operator 
can be moved inside the integrals to yield: 

) Jfe 2 d 4 Ax 2 M (A 4>(xi, yi, Ax)A^(x 2 ,y 2 , Ax)> W(xi, yi)W(x 2 , y 2 )dx!dy!dx 2 dy 2 (B.27) 

An expression for (A^(x ll yi, Ax)A<s{x 2 ,y 2 , Ax)) is therefore needed. Recalling the definition for 
A <t>: 

A^(x, y, Ax) = <t>(x — y) — 0(x + y) (B.28) 

The correlation function becomes: 


(A<£(xi,yi,Ax)A<£(x2,y2, 


Ax)) = ^(*1 - ^,yi)-^(xi + ^,yi) 






which can be expanded to be: 


(Ao(xi,yi, Ax)A<j>(x-i, jfc, Ax)) = (^ ^xi - ^ (*3 ~ 

~* (*1 “ 

-* (*i + ^r-w) * (**- x ,ys ) 

+4> (*i + ^ (*2 + IT’ 5 *)) 


(B.30) 


An expression for these correlation functions was found in Chapter III, Equation 3.32, to be: 


( <t>( x i<yi)<p{x2,y2)) 


_ _<L®2 ( 
" 2 V 


5/3 


V(»2 gi)^+(ya yi) 2 j +r*(o,o) (B.31) 


Therefore, by substitution of Equation B.31 and pulling out the constants, Equation B.30 becomes: 


(A<f>(x 1 , 3 / 1 , Ax) A*(x 2 , jfc, Ax)) = 


6.88 

2 r *' 3 


Ax Ax\ . , 2 

Z 2 + ___ Z 1 + _J 


+ /( Z2 ~ it ~ ri _ it) 

~^( Z2 _ it - Zl + ir) + ( y2 ~yi) 2 




"V (^ + x ~ 11 " t) + ( w-!/i ^ 


-f 


(B.32) 


Substituting Equation B.32 into Equation B.27, after some algebraic simplification, yields: 

(®?) = J[[[ [V(** ~ *t + Al ) 2 + (» - »0 J + \A* 2 ~ z 1 - ^*) 2 + (v* - yi) 2 


-2\/(x 2 - xi ) 2 + (1/2 -yi ) 2 


W(*i. yi)W r (x 2 , jft)dxidyidx 2 dy2(B.33) 
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Recalling Equation 2.42, the aperture function, Equation B.33 becomes: 


(<%) = k 2 d *±;2 r */3 M ^' 2 “ Xl + + “ VO * 1 + >/(*j-xi-A*) 2 + (ift-»i) a 

-2\/(x 2 — *l) J + (i/3 - yi) J 

xrect rect (~j^J rect rect dxidyidxjdyj (B.34) 

Equation B.34 can be written in dimensionless form by making the following substitutions: 

* = m - di ‘=&' 

= JJ 2 — d V. = |dV>< and 

Ax 

* = — (B.35) 

r o 

After some simplification, Equation B.34 then becomes 

/£»\ = 688 fd/rA 2 

\ •/ 2 4 2 5 / 3 ife 2 d 1 / 3 ro /3 V * / 

- 2 >/(&» ~ Si) 2 + (tk “ ^i) J j rect ^y^ rect ^y) rect ^y-^ rect ^y-^ 
xd^idV>id^ 2 d^ 2 (B.36) 

Note that all the terms inside the integrals are dimensionless, which allows solving Equation B.36 
to simply be an exercise in calculus. Therefore, the following variable substitution will be defined: 

and 5*=^y^- - ft = 5* - ^A* and ft = ££ + |a* 

At/> — ii) 2-^1 and ^ ^ = ^A^ and V'a = ^ Atf< (B.37) 

w 4 6 
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This choice allows the Jacobian to be one. Equation B.36 then becomes: 


Aa\_ 6-88 (d/r.y 

\ V 2*2 */W/3rJ /3 \ 3 ) 


M |( A4+2 5fe) ! + H , + |(^' 25 fe) , + a *1 , ' 2|a<,+A * , l* 

(5LJ*) „* t (H+i«) „cc (5fcJ*t) ^ (w+i*!) 


x d A£dAV>d££dE^ 


(B.38) 


Note that all terms that depend on Ei/> can be separate out. So, just for a moment, we can look 


only at the integral: 


J tM (2LzJ£ty ea (?*l il*t) ( n* 


(B.39) 


which is identical to Equation A.20 in Appendix A.l. Thus, from Equation A.23, this simply 
becomes: 

(B.40) 

Similarly, all terms that depend on can be separated from Equation B.38 to yield: 

out, 

Substituting Equations B.40 and B.41 back into Equation B.38 yields: 

Ap\_ 6.88 ( d/r o y 

• j 2 i 

x Jl UA£ + 2^-) +Ai> 2 \ +UA4-2—) + A4> 2 -2|A4 2 + AV. 2 | f 

x [tri [tri dA^dA^ (B.42) 
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Recalling the definition of the tri function, Equation A.27, setting the limits on the integration and 
performing some minor algebraic simplification, Equation B.42 becomes: 


<«>- 


6.88 




(nr) 7., f .?- |A{I)<2 - |A « (^ +2 3k ) !+ * 


K*- 2 ^) 2 


+ Atf 


-2|A£ 2 + AV> j |* 


dA4dA0(B.43) 


Which could be written as: 


where f v (s^f) is defined to be: 


^ = 0.13544/.(sfr) 


kW 3 rl /a 


(B.44) 


((Af + 2x) 2 + A xl> 2 * 


dA£dAV> (B.45) 


fv (x) = (x)~ 2 / / (2 - |A*|) (2 - \Arf>\) 

+ |((A£-2x) 2 + A^ 2 |* -2|(A£ 2 + Atf. 2 |* 

B.3 Derivation of Cross-Correlation of z-tili and s-tilt, (^0 t ■ 9,^ 

From Equation 4.23 


/a a \ _ l\~f*W x >yW{z,y)dzdy , ./y<£(x,y)W(z,y)dzdy 

V* ' t 'V ~ \ x M - + y - kP - 

\ 12 IT 

[ f J A<f(x,y, Az)W(z,y)dzdy ^ . /A^(x,y, Ay)W(x, y)dxdy j ^ ^ ^ 


A xkd 2 


+ y- 


A ykd? 


Performing the dot product yields: 


(St = Jf A<£(x,y,Az)W(z,y)dxdy^ ■ Jf x<£(z,y)W(x,y)dzdy^ 

+ (a yk<P Jf A<t, ( x 'y’ A y) W ( x ’y) dxd y'j ■ (Jffj* JJ !^(x,y)W(x,y)dxdy^(B.47) 
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Again, the first two integral pairs are essentially identical to the second two pairs. Hence, the 
second two integral pairs can be dropped and accounted for by multiplying the first two integral 
pairs by two. Combining the remaining integrals and pulling out the constants yields: 

£) = / JjH x 2^Hxi,yi,Ax)4>(x 2 ,y 2 )W(xi,yi)W(x 2 ,y 2 )dx i .dx 2 dyi<iy 2 ^ (B.48) 


Since <j> and Ao are the only random terms, the expected value operators can be moved inside the 
integrals to yield: 


(e t -e,) = 


A xk 2 d' 


x 2 (Atfi(xi, yi , Ax)<l>(x 2 , Vi)) W{x x , yi) W(x 2 , y^dxidx^ytdjft (B.49) 


Therefore, an expression for (A<p(xi, yi, Ax)<t>(x 2 , y 2 )) is needed. Again recalling Equation 2.13, 
this correlation function can be written as: 


{A<t>(xx,yi, Ax)<t>(x 2 ,y>)) = (j> ^ei - ^j-.yij <p{xi,yi) - <t> yijtf (* 2 , y 2 )^ (B.50) 


After using Equation 3.32, this simplifies to: 


{A<f>(xi,yi,Ax)4>(x 2 ,y2)) 


Axy , , 2 

X 2 -XX + — I +(y 2 -y i) - 




+ (y 2 -yi) (B.51) 


Substituting Equation B.51 into Equation B.49 yields: 


(*• *>) ][[[ * 2 W(x l ,y 1 )W(x 2 ,y 2 ) ^(*3-*! + ^) +(y 2 -yi) 2 

I -^-f ^ 

-y - Xi - + (y 2 - yi) 2 dz 1 dz 2 dyidy 2 (B.52) 
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Recalling the definition of the subaperture function, Equation 2.42, Equation B.52 becomes: 


' '*) " M ( 7 ) ( 7 ) ( 7 ) “« (f) 




+ (y2-yi) 2 ~ (*1 ~ x '~ *Y' S j +(v2-yx) 2 


xdxida:2dyidy2 


(B.53) 


As before, the integral in Equation B.53 will be normalized with the same variable substitution 
that was used in Section B.2, Equation B.35. After some simplification, Equation B.53 becomes: 

- g ) = M ( *' ct (I) ,ect (I) r “‘ (t) “* (f ) 

( \ 2 1 ^ 1 / \ 2 ^ 

* + + (^2-^i) 2 -K&—£ 1 -+(^2 _ ^i) 2 




(B.54) 


Again, following similar steps in as Section B.2, the variable substitution defined by Equation B.37 
is made, allowing Equation B.54 to become: 


(o,-o,) = 


12 • 6.88 d/r 0 

2525 / 3 ^^ 1 / 3 . 5/3 s 


m 


A* + ££ 


+ A^ 2 


x ^ + w-y + ^ -(«-£)' 


/E^ + iAV-\ JA _ 1A 


dA^d^dAtpd^ 


(B.55) 
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Examination of Equation B.55 shows that all terms that depend on E0 can easily be separated. 
Thus, just as in Section B.2, we can just look at the integral: 




(B.56) 


This is identical to Equation B.39 and simply becomes: 


/ (5*±J*t) „c. dE* = 2tri (f) (B.57) 


Hence, Equation B.55 becomes: 






2 i 

+ Ai/> 2 


xdA£dE£dA# 


(B.58) 


To aid in the visualization of the next steps, Equation B.58 will be written as: 


(‘ fl-sS*-** J> I-*)’— - («-*)•-’ 


xrect (5? rec t tri dA*dE£dAV< 

if«[ll 45+ 3^)’ + H* ■ 1^" + H* 

xrec t rect tri dA^dE^dA tp (B.59) 
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AS 


Figure B.l. Plot Showing Integration of an Odd Function over a Symmetric Period. 


In the first triple integral of Equation B.59, the terms depended on are easily separated to form 
the common result: 


(B.60) 


The E£ terms in the second triple integral can also be separated as: 

I (Hd&) ,ec« dE{ 


(B.61) 


While Equation B.61 is not as trivial as the others, examination of Figure B.l reveals that it is 
merely the integration of an odd function, E£, over a symmetric interval, and hence the integral is 
zero. 




(B.62) 


Therefore, the second triple integral in Equation B.59 goes to zero and, after substituting in Equa¬ 
tion B.60, Equation B.59 becomes: 
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(£•£)=; 


12-6.88 d/r 


2 3 2 5/3 jt 2 d i/3 r 5/3 s Jf 2 ^ (^d/fo) + | | (^ d/rj 

(¥)“(“) 


+ A^ 2 


xtn 


dA£dA^> 


(B.63) 


Recalling the definition of tri functions, Equation A.27, Equation B.63 becomes: 


<*■*)■ ^^ dJ r L L ^- |44l) < 2 " ** 

(«-wJ 


+ A t/» 2 


xdA£dAV> 


(B.64) 


Equation B.64 can be written as: 




kW 3 rl 13 


(B.65) 


where f c (s^) is defined as: 


fc (X) = z- 1 j 2 J A* (2 - |A£|) (2 - |A^|) 

x ||(A£ + xf + Att 2 1 ‘ - |(A4 - x) 2 + Atf 2 ’ j dA£dA^> 


(B.66) 
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